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I. INTRODUCTION 



In experiments, many hadrons are observed as resonances that decay via the strong in- 
teraction and have only a short life-time. On the theoretical side, the direct determination 
of the resonance parameters from QCD is afflicted with many difficulties since the compu- 
tation of resonance masses and decay widths is essentially a nonperturbative problem. An 
attractive way to extract the resonance parameters nonperturbatively from first principles is 
the use of lattice QCD. Among unstable hadrons, the case of the p-meson is ideal for lattice 
studies of a resonance for two reasons. First, in lattice calculations, the noise-to-signal ratio 
in the computation of a meson mass is proportional to e^ m ' M ~ nhT ' t , where vtlm is the meson 
mass under consideration, m n is the pion mass and t is a typical hadronic time scale. Since 
the p is one of the lightest mesons, the statistical error of its numerically computed mass 
can be well controlled. Second, the principle decay channel (with a branching rate of 99.9%) 
of the p-meson is to a pair of pions, which can be treated on the lattice very precisely. 

In the past, several lattice groups have undertaken efforts to study the p-meson decay. A 
first attempt was made to estimate the decay width from the p — > txtx transition amplitude []]- 
4|. This method relies on two assumptions: first, the energy gap between the ground and the 
first excited state (corresponding to p-meson and tttt states with the same quantum numbers) 
is small. Second, it is assumed that the hadron interaction is not large and the transition 
amplitude (p\irn) satisfies (p|vr7r) <C (p\p) l l 2 (ixix\ixix) l l 2 . An alternative method, which does 
not rely on these assumptions, is to extract the p-meson resonance parameters from the 
P-wave pion-pion scattering phase in the isospin 1 = 1 channel. The nonperturbative 
determination of the scattering phase is possible by using finite-size methods, which were 
originally proposed by Liischer in the center-of-mass frame (CMF) 5|-[9j and later extended to 



more general cases employing also a moving frame (MF) by Rummukainen and Gottlieb 101] 1 . 
Making use of these finite-size methods, two lattice studies fisl . [jjj] have been carried out 
to compute the p-meson resonance parameters. 2 These calculations mainly concentrated on 
the scattering phase at one or two energies for a single ensemble. In this way, however, the 
scattering phase can be extracted at only a small number of energies and it becomes difficult 



1 In Ref. [To| the determination of scattering phase is accomplished at the level of quantum mechanics. For 



the field theorectic derivations we refer to Refs 



11 



and, 



2 We note that recently another two lattice studies 15, lq were reported at Lattice 2010. 
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to map out the resonance region. 

In this work, we study the 1 = 1 pion-pion scattering system using three Lorentz frames: 
the CMF, the first MF with total momentum P = (2iv/L)e 3 (MF1) and the second MF with 
P = (27r/L)(ei + e 2 ) (MF2). Here, the denotes a unit vector in the spatial direction i 
and L is the spatial extent of the lattice. In each frame, we evaluate the P-wave scattering 
phase from the energy eigenvalues of the ground state and the first excited state. Using 
three frames allows us to obtain the scattering phase at six energies for each set of physical 
parameters considered without the need to go to larger lattices. Therefore, we think that 
our calculations have two advantages compared to the earlier works mentioned above. First, 
extracting the resonance parameters from six energies allows us to obtain more accurate 
results. Second, some of the scattering phases are calculated at energies that lie in the range 
[m p — r p /2,m p + Tp/2], allowing us to directly map out the resonance region. 

Our calculations are performed using the Nf = 2 maximally twisted mass fermion ensem- 



bles 



17 



19| from the European Twisted Mass Collaboration (ETMC) at a lattice spacing of 



a = 0.079 fm. The pion masses range from 290 up to 480 MeV, ensuring that the physical 
kinematics for the p-meson decay, m^jm 9 < 0.5, is satisfied. The computation of the p- 
meson resonance parameters at several values of the pion mass allows us to obtain the pion 
mass dependence of the resonance mass and decay width and hence to perform an extrapo- 
lation to the physical point. The benefit of using twisted mass fermions is that at maximal 
twist physical observables are automatically accurate to 0(a 2 ) in the lattice spacing, while 
the drawback is that isospin symmetry, although again an 0(a 2 ) effect for the observables 
considered in this work, is broken at nonzero values of the lattice spacing. As a result, for 
any value of a ^ the decay of p° to vr vr is allowed, while in the continuum limit isospin 
symmetry is restored and this decay is forbidden. In this paper we present a first calculation 
to extract the p-meson resonance parameters from three Lorentz frames and discuss the 
feasibility and accuracy achievable using this setup. Since here we use only one value of 
the lattice spacing, we cannot test for the possible effects of isospin breaking. We plan to 
come back to this issue in the future when we will analyze gauge field ensembles obtained 
at finer values of the lattice spacing. As it will turn out, we are not able to match the high 
experimental accuracy of the p-meson resonance parameters with our lattice calculation. 
Still, we consider this work an important conceptual study and the techniques used here will 
be useful for other resonances such as the A baryon. 
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II. METHOD 



A. Scattering phase 



In an elastic scattering system, the relativistic Breit-Wigner form (RBWF) for the scat- 
ring amplitude 
width T R is Q 



tering amplitude a\ with a resonance at a center-of-mass (CM) energy M R and with a decay 



-VsTr(s) _ p 2 



where Eqm is the CM energy and a>i is related to the scattering phase of the I th partial wave, 
5i, through ai = {e 2lSl — l)/2i. The RBWF corresponding to Si is then 

The p- resonance has quantum numbers I G (J PC ) = 1 + (1 ) and decays into two pions in 
the P-wave. A description of the scattering phase as a function of the Eqm is provided by 



the effective range formula (ERF) 21] 



^ = 9 h E ( J E 2 y P=y/&CM/i-mi, (2) 

which fits the experimental data well. In Eq. (J2]) 5\ is the P-wave pion-pion scattering 
phase, g pnn is the effective p — > mr coupling constant and m p is the p-meson mass. We 
remark already at this point that we will use the ERF also for our lattice calculations to 
fit the scattering phase, even when using pion masses that are larger than the physical one. 
Comparing Eqs. ([1]) and fl2]), we find that the ERF is a particular case of the RBWF if the 
parameters M R and r#(s) are chosen such that 

M R = m p , T R {s) = 9 -f^ P -. 

07T S 

The rho decay width V p can then be computed in the following way, 



r P = v R (s) 



9/mn Pp 

6tt m 2 p 



Pp = \ m l/ 4 ~ m l ■ ( 3 ) 



Thus Eqs. (T5]) and ([3]) allow us to extract m p and T p by studying the dependence of the 
pion-pion scattering phase 5± on Ecm- 
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B. Finite-size methods 



1. Center- of -mass frame 

A direct calculation of the phase shift from lattice QCD is possible by using a finite-size 
method established by Luscher js 9 ] . In this method, the phase shift is obtained from the 
energy eigenvalues of a two-pion system enclosed in a cubic box with spatial size L. 

In the CMF the possible energy eigenvalues for two noninteracting pions are given by 

E = 2^ml+p 2 , p=|p|, p=(2n/L)n, n G Z 3 . 
In the interacting case the energy eigenvalues are shifted, 



E = 2^ml +pi , p=(2n/L)q, 

where q is no longer constrained to originate from a quantized momentum mode. Because 
of the presence of the interaction, the energy eigenvalues deviate from those in the nonin- 
teracting case. It is exactly this deviation that contains the information of the underlying 
strong interaction and thus can be used to determine the scattering phase, as outlined next. 

In this paper, we concentrate on the energy eigenstates with energies E in the elastic 
region 2m 7T < E < Am^ with the two-pion system having the same quantum numbers as the 
p-meson. In the CMF these states transform as a vector (more specifically the irreducible 
representation V = T^) under the cubic group Oh- The corresponding finite-size formula 
connecting the energy E to the scattering phase 8\ is given by [9] 

3/2 

tanc^E) = J q , for T = Tf , (4) 

Am 1 ; q ) 

with the zeta function defined through 

^ooO;g 2 ) = —j= ^2 (l n l 2 _ <? 2 )~" • 



* n nez 3 
2. Moving frame 

Using a MF with total momentum P = (27r/L)d, d G Z 3 , the energy eigenvalues in the 
noninteracting case are given by 



E = J ml + p{ + J m 2 n + p 2 2 



where pi = \pi\ and denote the three- momenta of the pions, which satisfy the relations 

p, = (27r/L) ni , rij G Z 3 , pi + p2 = P . (5) 

In the MF, the center-of-mass is moving with a velocity of v = P/E. Using the standard 
Lorentz transformation with a boost factor 7 = 1/ a/1 — v 2 , the Ecm can be obtained as 

Ecu = T l E = 2^/ml+p* 2 , 

with CM momenta 

r = |p*i, p* = Pi = -p*2 = \r 1 (Pi-p2). (6) 

Here, we use the notation 

— 1 -1 P v 
7 P = 7 P|| + P± , P|| = — v > P± = P - Pll • 



v 



From inspecting Eqs. (EJ and ([6]) it can be seen that the p* are quantized to the values 

p* = (27r/L)n , n 6 — {n | n = 7 _1 (m + d/2) , for m G Z 3 } . (7) 

In the interacting case, the E CM is given by 

E CM = 2y/ml+p*i , p* = (2Ti/L)q. (8) 

From the energy shift between the noninteracting and the interacting situation, Eqm~ Ecm 
(or equivalently q 2 — |n| 2 ), one can compute the pion-pion scattering phase. 

In the MF1 (d = e 3 ), the energy eigenstates transform under the tetragonal group D^. 
The irreducible representations A2 and E~ are relevant for the pion-pion scattering states 
\iT7r,l = 1) in infinite volume with angular momentum 1 = 1. In this work, we calculate the 
energies associated with the A7[ sector. The formula converting the Ecm in a finite volume 



to the scattering phase in the infinite volume is given by Gottlieb and Rummukainen 10[ as 



tan 5i (Ecm) = —5 \- A , for T = , (9) 

Z d (l; g 2 ) + (2g- 2 /V5)2 2 d (l; g 2 ) 



with the modified zeta function 



^im\ s A ) ~ 2^ (\ n \ 2 -a 2 ' 

ngP d Vl 1 H ' 



and 

y lm (r) = r% m (Q r ) , y m (r) = r%- m (n r ) , 

where Q r represents the solid angle parameters (9, <fi) of r in spherical coordinates and the 
Y^ m are the usual spherical harmonic functions. 

In order to obtain more energies in the resonance region, we developed a second moving 
frame (MF2) with d = ei + The corresponding energy eigenstates transform under the 
orthorhombic group D^h- The irreducible representations B± , B^ and B% occur for the 
\ttit, 1 = 1) states in infinite volume. Here we focus on the B± sector. Our derivation of the 
corresponding finite-size formula for the MF2 results in 

tan 5i(E CM ) 



Z<x> - (q- 2 /V5)Zi + i(V3q-yVT0) {Z« 2 - ' 

for r = sr . (10) 



For more details, we refer the reader to Ref. 22 1. 



For brevity, we represent 2^(1; q 2 ) with the short notation Zf m in Eq. (flUI) . Using 
Eqs. (jl]), (Q and ffTOj) we can then convert a finite- volume determination of the Eqm into 
a calculation of the P-wave scattering phase Si. This is, of course, exactly the situation we 
are confronted with in a lattice calculation as performed here. 

C. Correlation matrix 

In the CMF, the value of the Ecu is directly given by the discrete energy eigenvalue 
E extracted from the large time behavior of the corresponding correlation function. In the 
MF, Ecm is related to E through the Lorentz transformation 

E 2 CM = E 2 - P 2 . (11) 

In order to calculate the energy eigenvalues E, we construct a 2 x 2 correlation function 
matrix through 



(**)(*) (^(0) ((tttt) (t) pt( )) 

^2x2{t) — X ' • (12) 

{pit) Ht(o)) {pit) pt (0 )> 
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1. tttt sector 



The tctc correlation function is constructed with the interpolating operators defined 
through 

(™)(*) = ^-J2MR) (n + (P/2 + Rp,t)7r-(P/2-Rp,t) 

ReG 

-7r"(P/2 + i?p,t)7r + (P/2- J Rp,t)) , (13) 

with the momenta on the lattice P and p taking discrete values 

P = (27r/L)d , p = P/2 + (27r/L)m , for d, m 6 Z 3 . 

Let us explain the notation we have used in Eq. ( 1T3"|) . The pion interpolating operator 
vr ± (q, t) is defined through 

*) = jj^ E e " tq ' X (^y A (x, f) , a = ±, , 

where r a denote the isospin Pauli matrices and ip the two-flavor quark fields. We also 
introduce the symmetry group G as the set of all lattice rotations and reflections R, under 
which the set of Pd defined by Eq. ([7]) is invariant 

G={A RneP dl VnGPd}. (14) 

In the CMF, MF1 and MF2, G is given by the cubic groups Oh, D^h and D 2 h, respectively. 
T is the irreducible representation of the group G, dr is the dimension of T and Xr{R) is 
the character of T. The average over all the operations R in the group G weighted by 
the coefficient Xr{R) projects out the scattering states that belong to the V representation. 
Finally, ]V G = 1. 

Given the momenta {P, p} and the representation T, one can construct the interpolating 
operators (7T7t) (t) using Eq. (|T3|) . Here we set T to be T±, and for the CMF, MF1 and 
MF2, respectively, so that the energy eigenstates |7r7r, T) in finite volume will approximate 
the P-wave scattering states \inr,l = 1) in infinite volume if one ignores states with higher 
angular momentum. In the CMF, the interpolating operator is given by 

(7T7T)(t) = 7T+ (y e 3,M ?r~ ( _ y e 3>M ~ n+ ("T^ 3 '*) ^ 
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In the two MFs, the operators are given in a unified form through 



(7T7T)(t) = 7T + (P, t) 7T" (0, t) - 7T + (0, t) 7l~ (P, t) , 

with P again the total threa-momentum of the scattering system. We can use these operators 
to measure the energy eigenvalues E from the corresponding correlation functions, convert 
E into Eqm by applying Eq. (jTTj) and then extract the P-wave scattering phase Si using the 
finite-size formulae listed above. 

2. p sector 

The interpolating operator for the neutral p-meson is constructed through a local vector 
current, 

p(t) = p°(P, f ) = J- £ e - P x (^(a • 7 ) y V) (x, t) , a • 7 = ^ a j7i , 

x ^ ' i=l 

where a indicates the polarization of the vector current. The direction of a is taken to 
be parallel to e3 in the CMF, e3 in the MF1 and ei + e2 in the MF2, respectively. This 
choice allows us to obtain a good signal-to-noise ratio for the off-diagonal matrix element 
(p(t) (7T7r)t(0)> in Eq. (02} . 

D. Extraction of energies 

By computing the matrix of correlation functions in Eq. ( 1T21) . we are able to isolate the 
ground state and first excited state in a clean way. This is of particular importance in 
the resonance region, where the avoided level crossing occurs and the first excited state is 
potentially close to the ground state. Such a situation renders the extraction of the ground 
state energy difficult when only a single exponential fit ansatz is used. Since we cannot 
predict a priori whether our energy levels will be close to the resonance region, we find it 
necessary to always use the correlation matrix to analyze our results. To extract the energy 
eigenstates, we follow the variational method [3] and construct a ratio of correlation function 
matrices as 

R(t, t R ) = C 2x2 (t)C^ 2 (t R ) , for t > t R , 
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Ensemble 


P 


afj, 


L/a 


77V 


m n /m p 


N 


A x 


3.90 


0.0085 


24 


480 


0.43 


176 


A 2 


3.90 


0.0064 


24 


420 


0.40 


278 


A 3 


3.90 


0.0040 


32 


330 


0.32 


124 


A A 


3.90 


0.0030 


32 


290 


0.30 


129 



TABLE I: Ensembles used in this work. We give the ensemble name Ai, the inverse bare coupling 
(3 = 6/<?q, the bare quark mass afi, the lattice size L/a and the value of in units of MeV. We 
also list the ratio m n /mp and the number iV of configurations used. 

where tn, the reference time slice, is assumed to be large enough such that the contributions 
to the matrix R(t,tn) from the excited states |n) with n > 2 can be ignored. 
The two eigenvalues R n (t,tu) (n = 1,2) of the matrix R(t,tn) behave as 

Rn(t, t R ) A n cosh (-E n (t - T/2)) , (15) 

where we assume that t is large enough (t > ^> 0) to neglect excited states but still far 
enough from the boundaries (t <C T/2) to ignore the unwanted thermal contributions as 



discussed in the case for the pion scattering length using twisted mass fermions in Ref. 23]. 



III. LATTICE CALCULATION 
A. Ensemble information 



The results presented here are from a sequence of ensembles with a lattice spacing of 
a = 0.079 fm. The pion masses range from = 480 to 290 MeV. At all pion masses the 
physical kinematics of m n /m p < 0.5 is satisfied, such that it is physically possible for the p- 
meson to decay into two pions. In our analysis we use two lattice sizes. The first corresponds 
to L — 1.9 fm with pion masses of = 480 MeV and = 420 MeV, i.e. ensembles Ai 
and A 2 in Table HI The second uses L = 2.5 fm with pion masses of m n = 330 MeV and 
m n = 290 MeV, i.e. ensembles A3 and A4 in Table D Additional information about the 
ensembles used is given in Table [land in Refs. fl7- 19|. 
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B. Sources 



To calculate the tttt correlation functions, we employ a stochastic method using Z4 noise 
sources £| (x) that are restricted to each three-dimensional time-slice with time t s . The 
sources Q (x) are also diluted in the color and spin indices, which are suppressed for sim- 
plicity. The index i runs from 1 to the number of stochastic noise sources. In this work 
we are able to achieve sufficient accuracy with just N s = 1 samples. Using the one-end 
trick 24|, we need to introduce two stochastic noise sources, e* q ' x £^(x) and Q b (x), for each 
factor of vr ± (q, t s ) in the correlation function. For the correlation functions in the MFs, we 
must account for two momentum modes (e iq ' x ^ s (x), q = and q = P). In the CMF there 
are three required momentum modes (e* q ' x ^ s (x), q = 0, q = (2ir/L)e 3 and q = — (27r/L)e 3 ). 
Since we place the source on all the time slices t s = 0, . . . , T — 1, we therefore perform T 
inversions for each configuration and each momentum mode. Note that the time extent of 
our lattices is chosen to be always twice the spatial extent. The correlator Cn(t) is then 
calculated through 

C n (t) = ((tttt) (f) (tttt)* (0)) = I £ ((tttt) (t + t s ) (tttt)* (t s )) . 

ts 

The rather large effort to generate propagators on all the time slices allows us to obtain the 
correlators with high precision, which is important to extract the desired energies reliably. 

In the calculation of the off-diagonal correlator, C*2i(t), the contraction of the quark 
fields leads to a three-point diagram. Since in this three-point diagram the two-pion fields 
are located at the same source time slice t s , we use the sequential propagator method to 
construct the correlator. We calculate C*2i(t) through 

c 21 (t) = wwt(o)) = ±£>(f + O0"r) t (*.)> , 

and again average the correlator over all time slices t s . For the second off-diagonal correlator 
Ci2(t), the two-pion fields are placed at the sink time slice t + t s , which would render the 
computation of C^it) more difficult. However, using the relation C^t) = C|iW) we § e ^ 
the off-diagonal matrix element C12 for free. 

For the p-correlator, C^t), we have performed a comparison between the stochastic 
source method and the point source method and found that the required computational 
effort to achieve a given signal-to-noise ratio is comparable. Historically, we started our 
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work with the calculation of the hadronic vacuum polarization tensor 25|. Since in that 
work we generated point source propagators for the ensembles listed in Table HI we just use 
the available propagators to construct the p-correlator 

c 22 (t) = (pHt + t s )p(t s )) , 

where now the source time slices, t s , are chosen randomly to reduce the autocorrelation 
between consecutive gauge field configurations. 

Because of the isospin symmetry breaking effects at nonzero lattice spacing in our maxi- 
mally twisted mass setup, the disconnected diagram for the neutral p-meson does not vanish. 
To address the disconnected contribution to the neutral p-meson, we need to generate, in 
principle, additional all-to-all propagators. However, the disconnected diagram correction 



has been studied in Ref. 



261 ] and has been found to be negligibly small, and hence we neglect 



it also here in the computation of the neutral p-correlator. For the same reason, we neglect 
the disconnected diagrams for the off-diagonal entries, where these contributions originate 
solely from the neutral p-operator. In the calculation of the correlator ((7r7r)(7T7r)'}, we are 
able to address these disconnected pieces, since we put stochastic sources on all the time 
slices. We find that the disconnected diagram makes an apparently small contribution to 
the correlator but adds a significant amount of noise, which would destroy the signal for 
the connected piece. Therefore we drop it from the tttt sector. To be clear, neglecting these 
disconnected contributions is not a genuine approximation but is simply ignoring lattice 
artifacts that would vanish in the continuum limit anyway. 



IV. RESULTS 



A. Energy eigenvalues 

In Fig. [Tjwe show our lattice results for R n {t,tji) {n = 1,2) in a logarithmic scale for 
the CMF, MF1 and MF2, as a function of time t together with a correlated fit to the 
asymptotic form given in Eq. ( 115)) . From these fits we then extract the energies that will 
be used to determine the scattering phase. Note that the slopes of \n(R n (t, £r)) are often 
very similar for n — 1 and n — 2, indicating that it is indeed essential to use the correlation 
function matrix. In order to extract the energies, we have to consider the two main sources 
of systematic error. One originates from the higher excited states and affects the correlator 
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FIG. 1: For the ensembles Ai (upper left), Ai (upper right), ^3 (lower left) and (lower right), 
the correlator R n (t,tji) (n = 1,2) as a function of t is shown. For each ensemble, from top to 
bottom the three plots present the lattice calculations in the CMF, MF1 and MF2, respectively. 
The solid lines are correlated fits to Eq. (I15|) . from which the energy eigenvalues E n are extracted. 
In each plot, the upper curve is n = 1 and the lower curve with the slightly steeper slope is n = 2. 

in the low-i region. The other arises from the unwanted thermal contributions that distort 
the correlator in the large-t region. By denning a fitting window [t m i n , t m ax] and varying 
the values of t min and t max , we are able to control these systematic effects. In practice, we 
set t min to be £r + 1 and increase the reference time slice tn to reduce the higher excited 
state contaminations. Besides this, we set t max to be sufficiently far away from the time 
slice t = T/2 in order that the fitting results are protected from the unwanted thermal 
contributions. The corresponding parameters t R , t min and t max used in this work are listed 
in Table IIVI All the ensembles shown in Fig. [T] visibly agree with the corresponding fit 
and lead to reasonable values of x 2 /dof. The x 2 /dof together with the fit results for E n 
(n = 1, 2) are also given in Table IIVI 
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B. Lattice discretization effects 



In the continuum limit, the Ecm is simply related to the energy spectrum E n through 
the Lorentz transformation of Eq. (II ip. However, on the lattice, the discretization effects 
explicitly break Lorentz symmetry and Eq. fllip is only valid up to discretization errors. 
Another discretization error arises from the continuum dispersion relation in Eq. flH]), which 
is particularly relevant for the finite-size methods used here. 



These two sources of systematic error have been studied in Ref. lOj, where the authors 
suggest to use the lattice modified relations 

cosh(£ CM ) = cosh(£ n ) - 2 sin 2 (P,/2) , n = 1, 2 , 

i 

cosh(£ CA/ /2) = 2sin 2 (p72) + cosh(m^) , p* = (2n/L)q , (16) 

instead of the continuum relations to reduce these discretization errors. Following this 
suggestion, we calculate the energy Eqm and the momentum p* from the energy eigenvalues 
E n using Eq. f fT6|) and then estimate the P-wave scattering phase 5i by employing p* in the 
finite-size formulae. The results for Ecm, P* and 8\ are given in Table IVl 



C. Extraction of resonance parameters 

From the Ecm we can now compute the P-wave scattering phases from six different 
energy levels, two from each of the three Lorentz frames employed. In order to extract the 
p-meson resonance parameters, we fit the results for the scattering phase to the effective 
range formula Eq. ($2§ and show the corresponding fits in Fig. [2J At the position where the 
scattering phase passes 7r/2, the resonance mass m p is determined. Additionally, the values 
of g p7 „T and hence T p are also evaluated from the fit. The corresponding results are given in 
Table HQ 

The finite-size methods are valid only for elastic scattering processes. In a situation with 
large enough energy, i.e. when Eqm > ^ m -K, the inelastic scattering channel will open and 
it is unclear how to the determine the scattering phase in such a case. Therefore, in our 
calculations we exclude results with energy E CM > 4:m n , which happened, fortunately, only 
for the excited state in the CMF for ensemble A±. 
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FIG. 2: We show for the ensembles A\ (upper left), A2 (upper right), A3 (lower left) and A4 
(lower right), the scattering phases calculated in the CMF, MF1 and MF2 together with the fits 
to the effective range formula Eq. ([2]). At the position where the scattering phase passes 7r/2, the 
resonance mass m p (denoted as ciMr in the graph) is determined. Through the fit, the coupling 
constant g pn7T and decay width T p are also extracted. 

D. Comparison with other results 

Using the resonance masses determined in the previous section, we show our values for 



m n . In order to compare 
as determined by the 



m p together with those of other groups in Fig. [3] as a function o 
these results, we scale m p and m n with the Sommer scale r 
groups individually. This avoids systematic effects when determining the lattice spacing 
from different observables and is most appropriate when one aims only at a comparison of 
results between different groups. We find a rather satisfactory agreement and attribute the 
mild variation among the groups with possible residual cutoff and finite-size effects in the 
various calculations, although a definite conclusion cannot be given here. 

We remark that our values of m p in physical units result from using the lattice spacing 
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m n (MeV) 


m„ (MeV) 


T p (MeV) 


9 pTTTV 


Ax 


480 


1118(14) 


39.5(8.2) 


6.46(40) 


A 2 


420 


1047(15) 


55(11) 


6.19(42) 


A 3 


330 


1033(31) 


123(43) 


6.31(87) 


A A 


290 


980(31) 


156(41) 


6.77(67) 



TABLE II: The results for the p-meson mass m p , the decay width T p and the effective p —> tttt 
coupling g pnn at pion masses ranging from 480 to 290 MeV. 

a=0.079 fm given earlier. This value of the lattice spacing was determined in Ref. by 
fixing the physical value of the pion decay constant f^. 



E. Quark mass dependence 

Having analyzed the ensembles listed in Table [J allows us to discuss now the quark 
mass dependence of the p-meson resonance parameters. There are several works using 
effective field theory (EFT) to describe the quark mass dependence of the p-meson resonance 



parameters 



34 



381 ] . The general structure of the pion mass dependence of m p and Y p can 



be written down as 

m p = M° p + C ml Ml + C m2 Ml + 0(M*) , 
T p = T° + C Tl Ml + C T2 Ml + 0(M*) . 

However, before using these formulae, it should be realized that m p and T p are not only 
statistically correlated, but also inherently related to each other, suggesting that the coef- 
ficients C m i and Cn (i=l,2) are not independent from each other. Therefore, in this work 



we will follow the strategy of Refs. 



39 



40] where m p and V p are considered as the real and 



imaginary part of the complex pole of the p-meson propagator. Hence, we introduce the 
complex pole parameter Z defined through 

Z = {m p - iT p /2) 2 . 

In this approach the power counting is given by the complex-mass renormalization scheme. 
Up to 0(q 4 ) in the chiral expansion, where q is a typical pion momentum, Z is given 
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l i ' 



• ETMC, n =2 

■ Graz, n f =2 
JLQCD, n =2 
PACS-CS, n=2+l 

► RBC-UKQCD, n =2+1 



0.5 



1.5 



(r n m ) 

v jt 



FIG. 3: p-meson mass as function of the pion mass squared, both scaled with r$. The p-meson 
resonance masses determined in our calculations (ETMC) are compared with those of the groups 



listed in the legend: chirally improved fermions (Graz) 
nonperturbatiyely improved Wilson fermions (PACS-CS) 



UKQCD) 



32!, 



28K overlap fermions (JLQCD) [23,1301, 



3l[ and domain wall fermions (RBC- 



331 ]. In order to be consistent, we include only the results of those groups for which 



we could readily find the values of ro evaluated at the same coupling and pion mass as is the 
p-meson mass. Also, note that only our calculation includes a proper treatment of the resonance 
nature of the p-meson. 



by [39|, |40] 



Z = Z x + C x Ml - 9 -^Z^Ml 



9, 



32tt 2 " V Ml 



Ml 



9 2 Ml 
16tt 2 Ml 



3 - 2 In 



Ml 
Ml 



2iir 



(17) 



where Z x = (M x — zr x /2) 2 is the pole of the p-meson propagator in the chiral limit, M 2 is 
the lowest-order expression of the chiral expansion for the squared pion mass and C x , g^p^ 
and g are coupling constants assuming real values. Using Eq. (]17|) to fit our results, we can 
determine the value of Z at the physical point, where it can be converted to the physical 
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resonance mass ?Tip, P hy and decay width r PiP hy. 

In practice, we perform the chiral extrapolation of Z in terms of the pion mass m % as 
extracted from the pseudoscalar correlation function as measured directly in the numerical 
calculations. By inserting the relation 

into Eq. (fT7|) . the expression for Z in terms of m n is given by 
Z-Z x + C x m„ ^if^Af 24^^^ 

^m^4-lV^4(3-2m4-2H ■ (18) 



19| with 



32?r 2 * V M 2 / 16tt 2 M 2 V M 2 
In Eq. f|T8|) the values of the input parameters F n and A 3 are taken from Ref. 

F n = i=/ = 86(1) MeV , In (A^/m^ phy ) = T 3 = 3.50(31) , 

where m^^hy is the physical pion mass. 

Before we perform a precise test of Eq. (ITS]) , we first confront our lattice results with a 
simplified fit ansatz to order 0(q 3 ), namely 

Z = Z x + C x ml- 9 ^Z^ml. (19) 

In the left panel of Fig. H] we plot the mass of the p-meson as a function of the square of 
the pion mass together with a fit to Eq. (fT9~]l . Using the fit to extrapolate to the physical 
point, our lattice result turns out to lie slightly high relative to the PDG value of the p-meson 
mass and shows a deviation of 1.9 a. 

In order to see whether higher-order corrections could reconcile our calculation with the 
experimentally determined p-meson mass, we also fit our lattice results to Eq. f lT8|) . All the 
fit results are listed in Table III II From the simplified fit to Eq. (|T9l) the lattice result of 



Af- 



g 2 /m /2Air is larger than the one suggested by EFT, which is g 2 pw /24ir = 3.4 GeV~ 
ter including the terms of 0(q ) in the fit, the uncertainty of the determination of ^ p7r /247r 
becomes, unfortunately, much larger and in fact g 2 p7r /247r cannot be determined in a statis- 
tically significant way. A similar situation happens in the determination of the parameter 
#7(16tt 2 M 2 ). The KSFR relation 



41 



42 
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Fit of Z to 


Eq. (JUD 


Eq. (USD 




0.821(24) 


0.850(35) 




0.171(31) 


0.166(49) 


M x 


0.756(24) 


0.803(47) 


r x 


0.190(35) 


0.179(58) 




6.42(45) 


4.9(2.0) 




9.8(1.5) 


10(12) 


ff 2 /(16vr 2 M x 2 ) 




0.01(1.09) 



TABLE III: The physical p- meson mass and decay width as extracted using Eq. (|19p and Eq. (|18[) . 
The values of m^phy, r Pi ph y , M x and T x are given in units of GeV and g 2 p7r /247r and <7 2 /(167r 2 M 2 ) 
are in units of GeV -2 . 

suggests that 5f 2 /(167r 2 M 2 ) takes the value of l/(327r 2 F 2 ) = 0.43 GeV" 2 . However, we are 
unable to determine g 2 / (167T 2 M 2 ) reliably from the fit. As can be inferred from Table ILTT1 
using a fit to Eq. (fl~8|) there is a 40% uncertainty in the determination of C x and a more than 
100% uncertainty in the determinations of both g 2 p7r /247r and g 2 / (167r 2 M 2 ). Proceeding with 
these results, we plot the mass of the p-meson as a function of the square of the pion mass 
together with the fit to Eq. ( jl~8l) in the right panel of Fig. HI At the physical point the 
result of m p is still high relative to the PDG value, suggesting that the pion masses used in 
the current calculations are too high for even the 0(g 4 ) extrapolations and that yet lighter 
quark masses will be necessary for quantitatively precise comparisons with experimental 
measurements of the p-meson mass. 

In Fig. El we plot the coupling g p7Tn as a function of the square of the pion mass and 
find that g p7T7T is practically independent of the pion mass. Moreover, the value of g p7T7T is 
consistent with the PDG value. This is not entirely unexpected. The coupling g pn7T , being 
dimensionless, is expected to be less sensitive to the pion masses and lattice spacings used 
in the calculation than the resonance mass m p is. In fact, whereas the accuracy of m p is 
currently systematically limited by the pion masses used in the calculation, the precision 
with which we can calculate g pn7T is clearly dominated by just the statistical errors of the 
current calculation. 
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m K 2 (GeV 2 ) m % 2 (GeV 2 ) 

FIG. 4: The />meson resonance mass as a function of the square of the pion mass. In the left 
panel, we fit the lattice results to Eq. (119p . In the right panel, we fit them to Eq. (|18|) . Note that 
these are combined fits to m p and T p (shown in Fig. [6]) . 

Eq. ([3]) shows that the decay width is determined from the fitted values of both m p 
and g p7T7T . Hence, we expect that it will reflect a combination of the aspects just discussed. 
In fact, in the chiral limit Eq. ([3]) reduces to T p = m p g 2 piT1T / '(487r). Thus the fact that 
m p overshoots the experimental measurement implies that T p will also be larger than the 
measured value. Additionally, the error of g p7r7r will be enhanced in T p leading to larger errors 
in the width than in the mass. These features can indeed be seen in Fig. El where we show 
the lattice results for V p as a function of the square of the pion mass together with the fit 
to Eq. (fT9l) in the left panel and with the fit to Eq. (1181) in the right panel. At the physical 
point, the decay widths are obtained as r P)P h y = 171(31) MeV using the fit to Eq. (|T9|) 
and as r PjP h y = 166(49) MeV using the fit to Eq. ffTB"]) . Both of the results are consistent 
with the PDG value F Pj pdg = 149.1(0.8) MeV within la. Note, however, that obviously the 
values determined from our lattice calculation are much less accurate than the one extracted 
from experimental measurements. Therefore, we consider the present work more as an initial 
study of how accurately resonance parameters can be extracted from nonperturbative lattice 
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calculations and not as a precise determination of these parameters. The results we have 
obtained here demonstrate that resonances can indeed be analyzed on finite lattices with 
numerical calculations. This is very promising, given the number of hadrons that appear in 
the physical QCD spectrum as resonances. 

10 
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FIG. 5: The effective coupling g p7Tn as a function of the square of the pion mass. 



V. CONCLUSION 

In this work, we have calculated the P-wave pion-pion scattering phase in the 1 = 1 chan- 
nel near the p-meson resonance region. We have performed our calculations at pion masses 
ranging from 480 to 290 MeV and at a lattice spacing of a = 0.079 fm. At all the pion 
masses, the physical kinematics for the p-meson decay, m 7T /m p < 0.5, is satisfied. Compared 
to previous calculations, we have pushed the techniques much farther forward by employ- 
ing three Lorentz frames simultaneously. This allowed us, in particular, to map out the 
energy region of the resonance without having to employ larger and more computationally 
demanding lattice calculations. 

Making use of Luscher's finite-size methods, we evaluated the scattering phase from six 
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FIG. 6: The p-meson decay width as a function of the square of the pion mass. The left panel 
shows the lattice results and the fit to Eq. (119j) . The right panel shows the fit to Eq. (I18p . Note 
that these are combined fits to T p and m p (shown in Fig. H]). 

energy eigenvalues per ensemble. In this way, we could fit the scattering phase with an 
effective range formula allowing us to extract the p- resonance mass m p , the decay width 
T p and the effective coupling g pnn . Taking the inherent relation between m p and V p into 
account, we have performed a fit to our results, obtained at four values of the pion mass, 
as a function of the complex parameter Z = (m p — iT p /2) 2 . This provided a means of 
extrapolation to the physical point. Even though our fit formulae are guided by EFT, our 
results are not precise enough to perform a thorough test of the fit ansatze. 

Keeping in mind the caveats just discussed, we quote for the p-meson mass m PiP hys = 
0.850(35) GeV and for the decay width r Pjphys = 0.166(49) GeV. When these values are 
compared to the corresponding experimentally measured quantities, it is clear that the 
lattice computations cannot yet match the experimental accuracy. Although a precise de- 
termination of resonance parameters on the lattice is still a challenge, our work serves as a 
next step in the attempt to understand the strong decays in a conceptually clean way. 
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Frame t R /a t min /a t max /a n x 2 /dof aE, 



CMF 7 8 

Ai MF1 9 10 

MF2 8 9 

CMF 8 9 

A 2 MF1 9 10 

MF2 9 10 

CMF 8 9 

A 3 MF1 10 11 

MF2 9 10 

CMF 8 9 

A A MF1 9 10 

MF2 9 10 
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2 1.90 0.4591(86) 

1 1.03 0.3363(14) 
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TABLE IV: Values of the energy eigenvalues for the ground state (n = 1) and the first excited 
state (n = 2) in the CMF, MF1 and MF2. In the table we list the ensemble number, the reference 
time tR, the lower and upper bound of the fitting window, t m ; n and i max , the fit quality x 2 /dof 
and the fit results for energy eigenvalues E n (n = 1,2). 
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Frame n aE n olEqm a P* <M°) 

1 0.4559(52) 0.1207(50) 137(3) 

CMF 

2 0.6584(90) 0.2686(57) 170(9) 
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CMF 
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A/TTTO 
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TABLE V: We give the P-wave scattering phase Si as extracted from the energies of the ground 
state and the first excited state in the CMF, MF1 and MF2. We list the ensemble number, the 
energies E n and Ecm, the momentum p* and the scattering phase 5\ (in units of degree). The 
single result marked by a star denotes that the corresponding Eqm is above the 4m n threshold. 
We therefore exclude that point from our calculations. 
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